Landau theory of superconducting phases in Th crystals 
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Abstract 

We consider the Landau theory of phase transitions for multiple superconducting phases in Th crystals. All 
possible phase transition sequences involving a single superconducting order parameter are found. These results 
may be applicable to PrOs4Sbi2. 
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Perhaps the greatest descriptive power and beauty 
of the Landau theory of phase transitions is revealed 
when successive phase transitions are described by a 
single, multi-component order parameter (OP). There 
are many well-known examples of this involving struc- 
tural OP's, especially among ferroelectric compounds. 
But although multi-component OP's are the hallmark 
of most unconventional superconductors, so far there 
are no confirmed examples of multiple superconduct- 
ing (SC) phases described by a single OP. PrOs4Sbi2, 
which displays two SC phases with different symme- 
tries, may be the first. 

The application of Landau theory to multiple SC 
phases described by a single OP was studied previ- 
ously [1] for various crystal symmetries, and applied to 
UPta . This material displays multiple superconducting 
phase transitions due to a two-dimensional SC order 
parameter, but the degeneracy of the order parameter 
is lifted by antiferromagnetic ordering. In PrOs4Sbi2, 
superconductivity occurs at Td = 1.85 K; an addi- 
tional phase transition is observed at Tc2 ~ 1.75 K 
[2,3] . Thermal conductivity measurements revealed the 
presence of nodes and a lowering of the symmetry of the 
gap function from four- fold to two- fold at Tc2 [4], which 
is highly suggestive of a symmetry lowering phase tran- 
sition involving the SC order parameter. In this arti- 
cle, we describe the phenomenological Landau theory 
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of multiple SC phases in Th crystals, which may be ap- 
plied to PrOs4Sbi2. 

Especially useful when the microscopic origin of su- 
perconductivity is unknown, the phenomenological ap- 
proach is based entirely on symmetry. Order parame- 
ters are classified according to representations of the 
crystallographic point group, in terms of which the 
Landau potential is expanded. In principle, all possible 
phases, their symmetries and phase diagrams, may be 
obtained in this way, as well as the positions and types 
of nodes of the corresponding SC gap functions. 

The point group Th has one-dimensional (ID), two- 
dimensional (2D) and three-dimensional (3D) repre- 
sentations, called A, E and T respectively, which may 
be even (for singlet pairing) or odd (triplet pairing). 
The phase diagram for the ID representation A in- 
volves only one SC phase, therefore we do not consider 
it any further. For the representations E and T, we re- 
strict our attention to phases connected by second or- 
der phase transitions, as suggested in at least one ex- 
periment [2]. Together with the restriction to a single 
OP, this sets strong constraints on allowed phase tran- 
sition sequences. 

The SC gap function is a 2 x 2 matrix, A(k) = 
iay'i/'(k) for singlet pairing and A(k) = i{d(^)B)ay for 
triplet pairing, where a = {ax,Sy,Sz) are Pauli ma- 
trices, t/'(k) is an even scalar function and d(k) is an 
odd pseudovector function. The gap in the quasiparti- 
cle energy spectrum in the singlet SC state is given by 
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A(k) = |^/)(k)|, while in the triplet state it is A±(k) = 
[|d(k)|2 ± |d(k) X d*(k)|]i/2. The functions V(k) and 
d(k) are expressed in terms of the components of the 
OP rji as 

V(k) = ^7yiVi(k), d(k) = ^7?idi(k), (1) 

i i 

where ipi(k) and di(k) are the basis functions for the 
even and odd representations of the point group, re- 
spectively [5]. 

The transformation properties of the OP are derived 
from the basis functions ipi and dj, which may be ar- 
bitrarily chosen to span the representation. For the 2D 
representation E, wo choose basis functions that are 
complex conjugate [6]. The Landau potential expanded 
up to sixth order in the 2D OP is 

F = ailmf + Imf) + + \V2f) + 2/32 1771^7721' 

+'n{vfvf + vhf) + I2i{rilr)f - vhf), (2) 

where a, fHi, /32, 7i, and 72 are phenomenological 
parameters. The Landau potential is invariant under 
K xU X G, where K is time reversal, U arc U{1) gauge 
transformations and G is the space group. Eq. (2) 
describes second order phase transitions between 
the normal state and the states (1,0) and {(f>i,4>2), 
as shown in Fig. la). The components of the state 
(01,02) have equal magnitudes and unfixed relative 
phases. The sixth order expansion is insufficient to 
describe the additional second order phase transition 
to the state (771,772), which has neither fixed relative 
magnitudes or phases. 

To obtain the third phase (771,772) it is necessary to 
consider additional higher order terms in the Landau 
potential. However, such an approach is rarely followed 
in practice. Instead, an effective Landau potential of 
an effective OP is used to analyse phase transitions be- 
tween the various SC phases whose symmetries satisfy 
group-subgroup relations. Second order phase transi- 
tions are indicated by the absence of third order terms 
of the effective OP. In our case, (771,772) has symme- 
try D2 , which is a subgroup of the symmetry group of 
(01, 02), D2 X K. Therefore the effective OP is ID and 
breaks time reversal symmetry, hence odd order terms 
in the effective Landau potential are prohibited. It fol- 
lows that the transition is second order. The symme- 
try of (1,0) is the group TiD-z) (sec Rcf. [6]), which is 
a supergroup of D2. The transition (1, 0) — > (771, 772) is 
described by an effective ID OP which is just 772. It is 
clear from Eq. (2) that there will be third order terms 
in the effective Landau potential of the form 772 ± 772^. 
Therefore the transition (1, 0) — > (771,772) is first order. 

For the 3D OP, real basis functions which transform 
under rotation as {x,y,z) are chosen for tpi and d, [6]. 
The free energy contains second and fourth order terms 
of the form 



a) Normal state 

(1,0) (cP„4>2) 

(171, 172) 

b) Normal state 

(1,1,1) (l,0,0)-^(|r]J,/ln2l,0) 

(l^il> \ni 0) ini, r)2, 0) (Inil, \ni) 

Fig. 1. Second order phase transition sequences for a) the 2D 
order parameter E and b) the 3D order parameter T. This 
diagram docs not show phases or transitions beyond two second 
order transitions from the normal state. Note that e — exp 
and {1, e, e^) has components with equal magnitude and phases 
fixed by s. 

h = |7?lP + \m? + 17?3P, n, h = |77ir + 1772^ + 1773^, 

^3 = rilrjf + ririf + rihf + c.c, 
and sixth order terms of the form 

Ilhh,hh,h = |7?ir+|7?2|« + |773|«,/5 = \m\\\ri2?- 

|7?3P) + |772[*(|7?3P - I771P) + |773|*(|77lP - W?), h,7 = 

[Vtvf + vhf + Vtvf ± (7?|77i*' + Vtrif + r^Uf)] + cc. 
These terms are sufficient to describe second order 
phase transitions between the normal state and the up- 
per four SC states shown in Fig. lb). As in the case 
of the 2D OP, higher order terms are needed to de- 
scribe the additional phase transitions shown in Fig. 
lb). Among the various SC phases satisfying group- 
subgroup relations, it is straightforward to determine 
which of them are connected by second order phase 
transitions, given their symmetry groups (see Ref. [6]). 
In all cases, the subgroup contains half the number of 
elements as the group, which means that the transi- 
tion can be described by a ID effective order parameter 
which has character —1 for half of the group elements. 
Third order terms in the effective potential are there- 
fore prohibited, and the transitions are second order. 

Properties of the gap function, especially nodes, may 
be determined by substituting the OP for each phase 
into Eq. (1). For PrOs4Sbi2, this analysis fails to ex- 
plain the four-fold to two-fold symmetry lowering of 
the gap function observed in Ref. [4] . This suggests that 
the transition at Tc2 is due either to a second SC OP, 
or else to an OP that is of a different origin altogether. 
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